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Abstract 

It is shown that in a class of maximal globally hyperbolic spacetimes admitting two local 
Killing vectors, the past (defined with respect to an appropriate time orientation) of any 
compact constant mean curvature hypersurface can be covered by a foliation of compact 
constant mean curvature hypersurfaces. Moreover, the mean curvature of the leaves of this 
foliation takes on arbitrarily negative values and so the initial singularity in these space- 
times is a crushing singularity The simplest examples occur when the spatial topology is 
that of a torus, with the standard global Killing vectors, but more exotic topologies are 
also covered. In the course of the proof it is shown that in this class of spacetimes a kind of 
positive mass theorem holds. The symmetry singles out a compact surface passing through 
any given point of spacetime and the Hawking mass of any such surface is non-negative. 
If the Hawking mass of any one of these surfaces is zero then the entire spacetime is flat. 



1. Introduction 

There are a number of general results in the literature on the properties of foliations 
by compact spacelike hypersurfaces of constant mean curvature (CMC hypersurfaces) in 
spacetimes which admit a compact Cauchy hypersurface. (See [16], [1] and references 
therein.) However, the only results which give criteria in terms of Cauchy data for the 
existence of such foliations covering more than a small neighbourhood of a given CMC 
hypersurface are restricted to special classes of spacetimes, all of which have high symmetry. 
The results of this paper also apply only to certain spacetimes with symmetry but represent 
a significant generalization, since they include for the first time spacetimes containing both 
matter and gravitational waves. The method used suggests that there is a close connection 
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between the question of global existence of CMC foliations and that of global existence of 
solutions of the Einstein-matter equations and such a connection helps to explain why it 
has up to now been necessary to make symmetry assumptions: we cannot understand the 
question of global existence of CMC foliations in a context where we do not understand 
the question of global existence for the Einstein-matter equations. When a CMC foliation 
exists in a spatially compact spacetime satisfying the strong energy condition it is unique, 
provided the exceptional case of flat spacetime is excluded. It provides an invariantly 
defined preferred time coordinate on spacetime. 

The spacetimes studied in the following are defined by two conditions. The first is 
that they be solutions of the Einstein equations coupled to certain matter fields and the 
second is that they admit a compact CMC Cauchy hypersurface which possesses a two- 
dimensional abelian group of local symmetries without fixed points. (The second condition, 
stated here informally, is made precise in the next section.) The simplest example of a 
symmetry of this kind is that where the compact Cauchy hypersurface is diffeomorphic to 
the torus T 3 = S 1 x S 1 x S 1 , with the symmetries given by the action of U(l) x U(l) 
acting on two of the three S 1 factors by rotations. As shown below, the mean curvature of 
the Cauchy hypersurface must be non-zero, except in the trivial case where the spacetime 
is flat. Without loss of generality, reversing the time orientation if necessary, it can be 
assumed to be negative. The main result of this paper (Theorems 5.1 and 6.1) is that if the 
spacetime is the maximal globally hyperbolic development of data with local U(l) x U(l) 
symmetry on a CMC Cauchy hypersurface then the entire past of this hypersurface can 
be covered by a CMC foliation, with the mean curvature taking all values in the interval 
(— oo,Hq], where H$ is the mean curvature of the initial hypersurface. In particular the 
initial singularity in these spacetimes is a crushing singularity in the sense of Eardley and 
Smarr[7]. The assumption made on the matter model is that it is either collisionless matter 
modelled by the Vlasov equation (Theorem 5.1) or a wave map with values in an arbitrary 
complete Riemannian manifold (Theorem 6.1). It is also shown that the CMC foliation 
can be extended so that the mean curvature takes on all values in the interval (— oo,0). 
Unfortunately this does not by itself suffice to show that the CMC foliation covers the 
entire future of the initial hypersurface. 

Special cases of this result are already known. The first is that of the Gowdy space- 
times on the torus. These are vacuum spacetimes with global U(l) x U(l) symmetry which 
satisfy the additional condition that the so-called twist constants vanish (see below). The 
result was proved in this case by Isenberg and Moncrief [12]. The second is that of solutions 
of the Einstein- Vlasov system with plane symmetry [17]. In the first of these cases there 
is no matter present, while in the second there are no gravitational waves. The results of 
this paper contain both these results as special cases. It should be noted that they go be- 
yond previous results even in the vacuum case in two ways: they require only local, rather 
than global symmetry and they allow non-vanishing twist constants. The essential new 
element in comparison with the cases considered previously is the occurrence of nonlinear 
hyperbolic equations which are coupled to the matter equations. These are treated with 
the help of methods introduced by Gu [11] in the study of wave maps and by Glassey and 
Strauss [10] in the study of the Vlasov-Maxwell system. 

As a by-product of this analysis, a theorem on the positivity of the Hawking mass in 
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spacetimes with local U(l) x U(l) symmetry is obtained. (Proposition 3.1). This says that 
the Hawking mass of any surface of symmetry in a spacetime of this type is non-negative 
and that if the Hawking mass vanishes for any one of these surfaces in a spacetime, then 
the spacetime is flat. 

2. Local U(l) x U(l) symmetry 

The spacetimes considered in the following are defined on manifolds of the form M = 
R x S, where S is a bundle over the circle S 1 whose fibre is a compact orientable surface 
F. Let p be the projection of the universal cover F onto F. Let g a p be a globally 
hyperbolic metric on M for which each submanifold {t} x 5 is a Cauchy hypersurface. S 
is covered by its pull-back to a bundle over R. Since R is contractible, the latter bundle is 
isomorphic to R x F. This is turn is covered in a natural way by R x F, which is simply 
connected. Hence the universal cover S of S can be identified with R x F and there is a 
natural fibre preserving projection corresponding to p. Let p be the associated projection 
of M = R 2 x F onto M. Define g a p to be the pull-back of g a p by p. Suppose that a 
two-dimensional Lie group G acts effectively on M by isometries of g a p in such a way 
that the orbits are the inverse images under p of the fibres of the bundle S. Each orbit 
with its induced metric is a simply connected Riemannian manifold of constant curvature 
and thus must be, up to a constant conformal rescaling, isometric to the standard metric 
on the sphere, the Euclidean plane or the hyperbolic plane. The isometry group of the 
sphere has no two-dimensional subgroups and thus the case F = S 2 is not possible. (If G 
is replaced by a three-dimensional Lie group then F = S 2 is possible and the spherically 
symmetric spacetimes studied in [17] are obtained.) The surfaces diffeomorphic to F which 
correspond to the fibres of S and whose inverse images are the group orbits will be referred 
to as surfaces of symmetry. These spacetimes will be said to have two-dimensional local 
symmetry. In the case where the orbits are isometric to the Euclidean plane, they will be 
said to have local U(l) x U(l) symmetry. 

Consider now any spacelike hypersurface S in the spacetime (M, g a p) which is a union 
of surfaces of symmetry. Choose one of these surfaces of symmetry and call it Fq. Let 7 
be an affinely parametrized geodesic of the induced metric on S which starts orthogonal to 
Fq. It is also orthogonal to all the other surfaces of symmetry which it meets. Taking all 
geodesies of this type and following them until they meet F again gives a smooth mapping 
(j) from F x / to S, where I is some interval. Let the parameter along the geodesies be 
chosen so that / = [0, 2%}. Let ifj denote the mapping which takes <f>(x, 0) to (f>(x, 2tt) and 
let ifi denote the lift of this mapping to a diffeomorphism between two inverse images of 
Fq in M defined by following geodesies in the covering space. The mapping ip maps the 
Killing vectors of Fq corresponding to the group action on one of the inverse images of Fq 
in M to those corresponding to another inverse image. These two inverse images can be 
identified with each other by an isometry, which is uniquely determined up to an element 
of the isometry group of the Euclidean or hyperbolic plane respectively by their induced 
metrics. In the case that Fq is isometric to the Euclidean plane i/jq is the composition 
of a linear mapping ifti with a translation ip2- The mapping if; must preserve the lattice 
of inverse images in F of a given point in F . These lattices are isometric, and so by 
using the freedom in identifying the two covering spaces, it can be assumed without loss of 
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generality that they are identified with each other. It follows that ipi can be represented 
as an element of GL(2, Z). When this element is not the identity it means in general that 
the topology of S is that of a non-trivial torus bundle over the circle. Let y A be periodic 
coordinates on Fo corresponding to Cartesian coordinates on Fo. Let (x, y A ) be Gauss 
coordinates based on F such that y A restrict to the previously chosen coordinates on F . 
Now let -B 4 = det(<7AB), where upper case Roman indices take the values 2 and 3. The 
metric takes the form: 

dx 2 + B 2 g AB (x)dy A dy B (2.1) 

where detg A B = 1- Let L be the length of a geodesic which starts normal to an orbit and 
ends when it intersects the same orbit again. Define a = 2ir/(f Q B~ 1 (x)dx) and: 

x' = af B- 1 (y)dy (2.2) 
Jo 

If x' is used as a coordinate and the primes omitted from the notation the metric takes 
the form: 

A 2 (dx 2 + a 2 g A Bdy A dy B ) (2.3) 

where A is a~ 1 B, reexpressed as a function of x' . The new coordinate x runs from to 27r. 
The functions <jab{x) satisfy the relation g A B {x + 2n) = n^n^gcDix), where are the 
components of a matrix in GL(2, Z). A satisfies A(x + 2n) = A(x). A similar analysis for 
the case that Fo is the hyperbolic plane would no doubt be more complicated and is not 
attempted here. However it apppears that, due to the fact that vector fields and tracefree 
symmetric rank two tensors on a surface of genus higher than one must have zeroes, in that 
case nothing will be obtained which goes beyond the spacetimes with hyperbolic symmetry 
already studied in [17]. 

The metric gAB can be parametrized in terms of two functions W and V in the 
following way: 

#22 = e^coshV, (733 = e _w/ coshV, (723 = sinhF (2-4) 

The values of V and W at x = and x = 2n are related by a diffeomorphism N which 
does not have a simple explicit form. There are several special cases which are of interest. 
Consider first the case where the matrix N with components is the identity. Then 
there is a natural action of {7(1) x U(l) on the spacetime and we have the case of (global) 
U(l) x i7(l) symmetry. A further specialization is given by the assumption that the 
reflections y A 1— > — y A are isometries of the spacetime metric for A = 2, 3. Spacetimes 
satisfying this condition will be called polarized {7(1) x {7(l)-symmetric spacetimes. They 
have the property that V = 0. The vacuum spacetimes of this class are the polarized 
Gowdy spacetimes [6]. The plane symmetric spacetimes studied in [17] have the property 
that W = V = 0. When is not the identity, there are two qualitatively different cases. 
If N is diagonalizable and not the identity, then either it is minus the identity, or the 
two eigenvalues are distinct. If it is minus the identity then S has a two- fold cover which 
is a torus and for our purposes is essentially the same as when N is the identity. When 
the eigenvalues are distinct the manifold S admits a geometric structure of type Sol in 
the sense of Thurston [18]. The metrics obtained in that case include ones which are of 
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Bianchi type VI . There is a polarized case, where reflections in the eigendirections of N 
are supposed to be isometries of the metric g a p on the universal covering space. If N has a 
non-standard Jordan form then, by passing to a two-fold cover if necessary, we can assume 
that these eigenvalues are equal to unity. The resulting manifold S admits a geometric 
structure of type Nil [18]. The metrics obtained in that case include those of Bianchi type 
II. 

Lemma 2.1 Let (M,g) be a non-flat spacetime with local U(l) x U(l) symmetry having a 
symmetric constant mean curvature Cauchy hypersurface and satisfying the dominant and 
strong energy conditions. Then given any point p on the Cauchy surface there exists an 
open neighbourhood U of p and a smooth local diffeomorphism of I x [0, 2tt] x T 2 onto 
U for some interval I such that: 

(i) if 4>(t, x±, yi) = <f)(t, X2, 2/2) then x\ = and X2 = 2tv or vice versa 

(ii) for each t e I the set 4>({t} x [0, 2tv] x T 2 ) is a hypersurface of constant mean curvature 
t. 

(Hi) the pull-back of the metric under 4> has the form 

-a 2 dt 2 + A 2 [(dx + (5 1 dt) 2 + a 2 g A B(dy A + [5 A dt)(dy B + (5 B dt)} (2.5) 

The functions a, (3 a , A and qab depend on t and x and c/ab has unit determinant. They 
satisfy a(t, 2tt) = a(t,0), A(t,2n) = A(t,Q), ^(t, 2tt) = /^(t.O) = 0, (3 B (t,0) = 0, 
gAB{t, 2tv) = n^n^c/cDit, 0), where is an element ofGL(2, Z). The quantity a depends 
only on t. 

Proof It is a standard fact that, in a non-flat spacetime satisfying the strong energy con- 
dition, a neighbourhood of a compact CMC hypersurface with non-zero mean curvature 
can be foliated by constant mean curvature hypersurfaces and that the mean curvature 
of these hypersurfaces can be used as a time coordinate. If the U(l) x £/(l)-symmetry 
is global then it follows from the uniqueness of CMC hypersurfaces that they are unions 
of surfaces of symmetry. If the symmetry of the data is only local then some more care 
is needed, but since an almost identical argument has been given in [17] the details are 
omitted here. The mean curvature of the Cauchy hypersurface in the assumption of the 
lemma cannot be zero. To see this consider the Hamiltonian constraint 

R - k ab k ab + (trk) 2 = 16np (2.6) 

When the mean curvature is zero this implies that the scalar curvature R is non-negative. 
The topology of the Cauchy hypersurface is such that any metric with non-negative scalar 
curvature must be flat. For its universal cover is diffeomorphic to R 3 . This implies 
([14], p. 324) that the Cauchy hypersurface admits no metric of positive scalar curvature 
and it is well-known that a compact 3-manifold satisfying the latter condition admits 
no non-flat metrics of non-negative scalar curvature. Hence the induced metric on the 
Cauchy hypersurface is flat and, from the Hamiltonian constraint the second fundamental 
form and the energy density are zero. It follows from the dominant energy condition 
that the spacetime is vacuum everywhere and uniqueness in the Cauchy problem for the 
vacuum Einstein equations shows that the spacetime is flat. Since the spacetime is by 



5 



hypothesis non-flat, it follows that a maximal hypersurface is impossible. Let U be an 
open neighbourhood of the Cauchy hypersurface covered by a CMC foliation and let t be 
the function on U which is equal to the mean curvature of the leaf of the foliation on 
which the point lies. Choose a surface of symmetry Fq in the initial hypersurface t =const. 
and identify this with surfaces in the other hypersurfaces t =const. by means of geodesies 
which start on Fq orthogonal to the Cauchy hypersurface. Construct a mapping on each 
hypersurface t =const. from [0, 2tt] x T 2 in the way described above. Putting together 
these mappings for all values of t occurring in the foliation of U gives the mapping whose 
existence is asserted by the lemma. □ 

3. Estimates for the Hawking mass and area radius 

In this section certain general estimates for spacetimes with local U(l) x U(l) symme- 
try are derived. A solution of the Einstein constraint equations consists of a 3-dimensional 
manifold S and a Riemannian metric h ab , a symmetric tensor k ab , a real- valued function p 
and a covector j a on M which satisfy the Hamiltonian constraint (2.6) and the momentum 
constraint 

V a k ab -V b (trk)=8ivj b (3.1) 

In this paper it is always assumed that the dominant energy condition holds and this implies 
that p > \j a \. Suppose now that S is covered by a Gaussian foliation. In other words, if Fq 
is a fixed leaf of the foliation, any other leaf is obtained by going a fixed distance along the 
geodesies which start normal to Fq. If we think of S as being embedded in spacetime, then 
the resulting embedding of each leaf F in this spacetime defines various geometrical objects 
on F, as is always the case for an embedding of pseudo-Riemannian manifolds. We present 
the definition of these objects in terms of such an embedding, but in fact they are uniquely 
defined by the initial data. There is a preferred orthonormal basis of the normal bundle 
of F in spacetime, where the first vector is normal to S and the second vector tangent 
to S. These vectors are defined uniquely up to sign by this condition. The geometric 
objects defined on F by the embedding are then the induced metric, a second fundamental 
form associated to each normal vector and a 1-form, which is the representation of the 
normal connection in the given normal basis. The two second fundamental forms will be 
denoted by kab and Xab respectively and the 1-form representing the normal connection 
will be denoted by i]a- (Here upper case Roman indices are used for objects intrinsic to 
F. Indices of objects of this kind are raised and lowered using the induced metric qab 
and its inverse.) In a vacuum spacetime the freedom in t\a consists of just two spacetime 
constants. These are the twist constants referred to in the introduction, whose vanishing 
is one of the defining conditions of Gowdy spacetimes. If u is an arc length parameter 
along the normal geodesies, the constraints can be written in the following form: 

d u {tx\ + tr/c) = H(tr\ + tm) + V A r] A + K 

- 8ir(p + J) - |(trA + tr«) 2 - \{\ AB + k AB )(~X AB + h AB ) - ??V (3.2) 
d M (trA - tr«) = -H(tr\ - tr«) - V A T] A + K 

_ 87r (p - J) - f (trA - tr«) 2 - \{\ AB - k AB ){\ AB - ~kab) - ??V (3.3) 

d u rj A = -(tr\)r] A - V B k A b - 8jvj A (3.4) 
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Here H is the trace of the second fundamental form k a b (i.e. the mean curvature), K 
is the Gaussian curvature of F, Rab and Xab are the trace-free parts of kab and Xab 
respectively and J is the contraction of the unit normal vector to F in S with j a . This way 
of writing the constraints generalizes an approach used by Malec and Murchadha[15], 
for spherically symmetric asymptotically flat spacetimes, by the author [17] for spatially 
compact surface symmetric spacetimes and by Chrusciel[5] for vacuum spacetimes with 
{7(1) x i7(l) symmetry. In the following these equations are only used in the case of local 
U (1) x {7(1) symmetry. It should, however, be noted that the form of the equations suggests 
that there may exist an analogue in the general case. The terms which cause difficulties 
in general are those containing derivatives tangential to the foliation by surfaces F. These 
are of the form V a t]a and K. The integral of the first of these over F is zero, while the 
integral of the second is a constant only depending on the topology as a consequence of 
the Gauss-Bonnet theorem. Thus integrating over F eliminates the tangential derivatives 
from equations (3.2) and (3.3). 

Consider now the case of local {7(1) x {7(1) symmetry, with the Gaussian foliation 
being that by surfaces of symmetry. Let 9 = trA — tin, 9' = trA + tr«. These are the 
expansions of the two families of null geodesies orthogonal to F. Let the area radius r be 
the square root of the area of F. The Hawking mass is defined by m = — ^rV a rV Q r. In 
this case V a t]a = K = and equations (3.2) and (3.3) become: 

d u 9 = -H9 - P (3.5) 
d u 9' = H9' - P' (3.6) 

where the quantities P and P' are non-negative. It is also useful, following [15], to write 
these equations in the alternative form: 

8 u {r9) = -Q- ~^(9 2 r 2 + A9Hr 2 + 9r(9r - 9'r)) 

\ (3-7) 

d u (r9') = -Q' - — {9 2 r 2 - 49' Hr 2 + 9'r (9'r - 9r)) 
Ar 

where the quantities Q and Q' are non-negative. Consider now a symmetric Cauchy 
hypersurface S, i.e. one which is a union of surfaces of symmetry. Denote the maximum 
value attained by r9 and r$' on this hypersurface by M_|_ and the minimum by M_. Let 
xq be a point where the maximum is attained and suppose without loss of generality that 
9(xq) > 9'{xq). Since xq is a critical point of r6>, it follows from (3.7) that at that point 
either r9 < or 

9 2 r 2 + (4Hr)(9r) < (3.8) 

It follows that M_|_ < 4| Hr | . Similarly, M_ > —A\Hr\. These inequalities show that 9 and 
9' can be bounded in modulus by 4|i7|. The Hawking mass is related to the area radius and 
the expansions by — 2m/r = \r 2 99' . Thus in a spacetime with local U (1) x {7(1) symmetry 
which is foliated by compact CMC hypersurfaces with the mean curvature varying in a 
finite interval (£1,^2) and which satisfies the dominant energy condition, if r is bounded 
then 2m /r is bounded. These equations can also be used to prove a kind of positive mass 
theorem. 
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Proposition 3.1 Let (M,g) be a spatially compact spacetime with local U(l) x U(l) sym- 
metry which satisfies the dominant energy condition. Then the Hawking mass of each 
surface of symmetry is non-negative and if the Hawking mass of any surface of symmetry 
is zero the spacetime is flat. 

Proof The proof is similar to that of Lemma 2.4 of [17]. If m vanishes on some surface 
F then 9 or $' is zero there. Suppose without loss of generality that it is 9. Then it can 
be concluded as in the proof of Lemma 2.4 of [17] that 9 and P vanish on any symmetric 
compact Cauchy hypersurface containing F. When 9 is zero the other expansion 9' is given 
by the rate of change of r along the compact Cauchy hypersurface. Hence 9' must vanish 
somewhere and, repeating the previous argument, 9' and P' vanish on the whole Cauchy 
hypersurface. Looking at the explicit forms of P and P' shows that p = 0, kab = 0, 
^ab = and t\a = on the Cauchy hypersurface. The vanishing of \ab implies that c/ab 
is independent of x. It follows that a linear transformation with constant coefficients of 
the coordinates y A can be used to reduce the metric cjab on a given Cauchy hypersurface 
to the form Sab- This, together with the vanishing of kab and tja, shows that the initial 
data are plane symmetric. The subset of spacetime where m = is closed. Because 
of the possibility of deforming spacelike hypersurfaces, it is also open and must be the 
whole spacetime. Hence the spacetime is plane symmetric and applying Lemma 2.4 of 
[17] shows that it is flat. If the spacetime is not flat then it follows that 9 and 9' can 
never vanish. If they had opposite signs then this would mean that the gradient of r was 
everywhere spacelike and hence that the restriction of r to a Cauchy hypersurface was 
strictly monotonic. This is clearly impossible, since this restriction must have a critical 
point somewhere. Hence 9 and 9' have opposite signs, the gradient of r is timelike and the 
Hawking mass is positive. □ 

The timelike vector V a r is past-pointing. For otherwise 9 would be negative and 9' positive. 
Integrating (3.5) from to 2n on a hypersurface of constant time would then imply that H 
was positive somewhere, contrary to what has already been assumed. It follows that r is 
non-decreasing to the future along any timelike curve and that its value at any point with 
time coordinate t\ is bounded by above by its value on the hypersurface t = ti if t\ <ti- 

Let n a denote the unit normal to the surfaces of symmetry in the hypersurfaces 
t =const. and define K = k a bn a n b . Then, with respect to the coordinates introduced 
in Lemma 2.1, some of the field equations take the following explicit forms: 

dl(A 1/2 ) = -\A^ 2 \l(K - \t) 2 - It 2 + 2 VA v A + k AB kab + \ AB \ab + 167rp] (3.9) 



d 2 x a + A~ 1 d x Ad x a = aA 2 [^(K - \t) 2 + \t 2 

+ 2 V av A + k AB k ab + 4n(p + tiS)] - A 2 (3.10) 

d x K + 3A~ 1 d x AK - A~ 1 d x At - k ab X A b = 8nJA (3.11) 

d x ^ = -a^dta + \a(?>K - t) (3.12) 

d t a = al-d^ 1 + \a{ZK - t)\ (3.13) 

d t A = -aKA + d^A) (3.14) 



These equations have been written in a form which makes as clear as possible how they 
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differ from the form they take in the special case of plane symmetric spacetimes. The 
differences are not very great and in particular equations (3. 12)- (3. 14) are identical to the 
corresponding equations (2.6)-(2.8) in [17]. In terms of these variables the expansions are 
given by 9 = 2A~ 2 d x A -t + K,6' = 2A~ 2 d x A + t-K. 

On an interval where H is bounded the quantities 9 and 9' and K = t + \{9 — 9') 
are bounded. On the other hand, It follows from the lapse equation (3.10) that a < 3/t 2 . 
Integrating equation (3.13) in space shows that on a finite time interval a and a -1 are 
bounded. Putting this back into the integrated equation shows that dta is bounded. 
Equation (3.13) then implies that d x f3 x is bounded. Equation (3.14) can be rewritten as 

t (log4) -/^(logA) = -aK + dx/3 1 (3.15) 

Together with the bounds which have just been derived this implies that A and its inverse 
are bounded. Since r = aA it follows immediately that r and its inverse are bounded. The 
inequalities just obtained serve as a replacement for the bound for m _1 obtained at the 
corresponding point in the argument in [17]. The argument of [17] would apparently not 
work in the present case because tja makes a contribution to the equation for V a w which 
has the wrong sign. The argument used here also has the advantage that it only requires 
the matter to satisfy the dominant and strong energy conditions and no assumption on 
the positivity of the pressure is needed. For this reason it applies to more general matter 
models and in particular to situations where an electromagnetic field is present. 

The following theorem can now be proved: 
Theorem 3.1 Let a solution of the Einstein equations with local U(l) x U(l) symmetry 
be given and suppose that when coordinates are chosen which cast the metric into the form 
(2.5) with constant mean curvature time slices the time coordinate takes all values in the 
finite interval (£1,^2) • Suppose further that: 

i) the dominant and strong energy conditions hold 

ii) t 2 < 

Then the following quantities are bounded on the interval (ti,t 2 ): 

a, d x a, A, A -1 , d x A, K, 1 , a, a -1 , d t a (3.16) 
d t A, d x (3 x (3.17) 

Prooflt has already been shown that a, A, A -1 , a, a -1 , K and dx/3 1 are bounded. The fact 
that 9 and 9' are bounded implies that A' is bounded. The boundedness of d x fi x and the 
fact that P 1 vanishes at one point show that fl 1 is bounded. Equation (3.14) gives a bound 
for dtA. Integrating equation (3.9) over the circle and using the bounds obtained already 
shows that f Q n p and f Q n {2r\AT\ A + k ab kab + ^ AB ^ab) are bounded. By the dominant 
energy condition it follows that j and trS are bounded. Finally, integrating (3.10) 
starting at a point where d x a = gives a bound for d x a. □ 
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4. Estimates for the hyperbolic and Vlasov equations 

The field equations which are used to control W and V are hyperbolic. These quan- 
tities may be thought of as describing gravitational waves. The fact that these equations 
are coupled with the matter equations and are themselves nonlinear means intuitively that 
the waves interact with the matter and with each other. The equations will be written 
in terms of a 2+2 split of the metric. Here lower case Roman indices refer to objects 
which live on the quotient of spacetime by the symmetry group. Indices of objects of this 
kind are raised and lowered using the metric g a b on the quotient space and its inverse. The 
equations are: 

V a (r 2 V a W) = -r 2 tanhVV a WV a V - r 2 (coshV)- 1 [e- w T 2 2 - e w T 33 

-\{e- w {mf-e w {m) 2 )] (4-1) 

V a (r 2 V a V) = r 2 cosh V smhVV a WV a W - 2r 2 (cosh ^) _1 [( T 23 - \~h AB T AB g 23 ) 

-Uv2V3-U~ hAB VAVB)~92s)] (4.2) 

The derivation of these equations is lengthy and it proved useful for this purpose to make 
use of the calculations of Kundu [13]. Let Sw an d Sy denote the right hand sides of 
equations (4.1) and (4.2) respectively. It will now be shown that the modulus of each 
of these quantities can be bounded by a constant multiple of the expression p + t)aV A + 
kabk ab + Xab^ AB • It then follows from what was said in the last section that under 
the hypotheses of Theorem 3.1 the L 1 norms of Sw and Sy in space are bounded by a 
constant which does not depend on time. For this purpose it is necessary to calculate Xab 
and Rab explicitly in terms of W and V. 

\ AB \ AB = iA- 2 (cosh 2 VW 2 + V 2 ) (4.3) 
k AB k AB = i«- 2 [cosh 2 V(W t - pW x ) 2 + (V t - pV x ) 2 \ (4.4) 

This shows that the first term on the right hand side of each of the equations (4.1) and 
(4.2) can be bounded by Xab^ AB + kabk ab . To bound the other terms on the right hand 
side of (4.1) and (4.2), define an orthonormal frame on each orbit by: 

e 2 = (Aa)-\e- w/2 cosh(V / 2) 8 /Oy 2 - e w/2 sinh(V/ 2)0 /8y 3 ) (4.5) 

e 3 = (Aa)-\-e- w / 2 smh(V/2)d/dy 2 + e w l 2 cosh(y/2)a/a ? / 3 ) (4.6) 

Then 

e~ w/2 d/dy 2 = Aa[cosh(y/2)e 2 + sinh(y/2)e 3 ] (4.7) 

e w/2 d/dy 3 = Aa[sinh(y/2)e 2 + cosh(y/2)e 3 ] (4.8) 

The components of the covector t\a expressed in an orthormal frame can be bounded in 
terms of i] A i]A- Thus if the latter expression is bounded it follows that the components of tja 
expressed with respect to the basis (e~ w / 2 d/dy 2 , e w / 2 d/dy s ) can be bounded by a constant 
multiple of cosh(V/2) or, equivalently, by a constant multiple of (coshF) 1 / 2 . This means 
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that e~ w l 2 r}2 and e w / 2 i] 3 can be bounded by an expression of the form CV^^cosh V) 1 / 2 
for some constant C. This allows the expressions on the right hand side of equations (4.1) 
and (4.2) containing tja to be bounded in modulus by a constant multiple of tjaT) A . The 
terms involving the energy-momentum tensor can be handled in a very similar way. The 
dominant energy condition implies that the components of the energy-momentum tensor 
in an orthonormal frame are bounded in modulus by p and using (4.7) and (4.8) allows 
this to be translated into a bound on the matter terms on the right hand side of equations 
(4.1) and (4.2) in terms of p. 

Lemma 4.1 Under the hypotheses of Theorem 3.1 the quantities W , V , tja, (3 a , d x (3 A are 
bounded. 

Proof Choose some £3 in the interval (ti,t 2 ) and let (£4,0:4) be some point of the quotient 
manifold M parametrized by t and x with £4 < £3. (The case £4 > £3 is similar.) The 
equations (4.1) and (4.2) have the same characteristics. These are the null curves of the 
metric defined on M. Let 71 and 72 be the two characteristics passing through (£4, X4) and 
let (£3,2:5) and (£3,^) be the coordinates of the points where they meet the hypersurface 
£ = £3. The left hand side of equation (4.1) has the form of a divergence. Applying Stokes' 
theorem to the triangular region T bounded by 71, 72 and the curve £ = £3 gives the 
identity: 

,t 3 r t 3 



S w aAdtdx= I ' (r 2 DW/Dt)(t,<y l (t))dt+ f * (r 2 DW/Dt)(t,<y 2 (t))dt 

J td J td_ 

' r 2 (W t -[3 l W x )(t 3 ,x)Adx 



X 5 



and hence, after integration by parts: 

(r 2 W)(t 4 ,x 4 ) = i[(r 2 W0(£ 3 ,x 5 ) + (r 2 W)(t 3 , x 6 )} 

- \ f 3 (2rWDr/Dt)(t, ll (t))-^ f * (2rWDr/Dt)(t, l2 (t)) 

J td_ J td_ 



(4.9) 



- \ I \ 2 {W t - (3 1 W x )(tz,x)Adx-\ ( S w aAdtdx 

Jx 5 JT 



Here D/Dt denotes a derivative in the direction of the characteristic along which the 
integration is carried out, with this characteristic being parametrized with respect to £. In 
other words, for any function /, Df /Dt = d/dt(f{^{i))). Most of the quantities in (4.9) 
are already known to be bounded. This is in particular true of Dr/Dt. Thus the following 
inequality holds: 



|r 2 W(£ 3 -£)|U<C 




/•72(*3-s) \ 

r 2 W(£ 3 -s)|U+ / \S w (s,x)\dx)ds 



(4.10) 



Since ^ \Sw(t, x)\dx is known to be bounded and under the given hypotheses the number 
of times the characteristics can go around the circle between £ = £1 and £ = £3 is bounded, 
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it follows from Gronwall's inequality that W is bounded. The same kind of argument shows 
that V is bounded. The remaining conclusions of the lemma are simple consequences of 
the boundedness of W and V, as will now be shown. The momentum constraint implies 
that: 

d x (A 2 VA) = 87iA 2 j A (4.11) 

which means that (A 2 rfA)(t, x\) — (A 2 r]A)(t, X2) is bounded independently of t, x\ and 
X2- On the other hand the boundedness of r] A r]A(x)dx together with that of V and 

W shows that \A 2 t]a(x) \dx is bounded. These two facts together show that t\a is 
bounded. The definition of the second fundamental form gives the equation: 

9ABd x (3 B = 2aA- 1 a- 2 VA (4.12) 

This means that d x /3 A is bounded and, remembering that by definition f3 A (0) = 0, this 
implies that (3 A is bounded. This completes the proof. □ 

Everything which has been done up to now consists of obtaining bounds for parts 
of the geometry using nothing about the matter model except the dominant energy and 
non-negative pressures conditions. Now the special case of the Vlasov equation will be 
considered. In this class of spacetimes the Vlasov equation for particles of unit mass takes 
the following form: 

df/dt + (aA-^v^v ) - P^df/dx + F^f/dv 1 = (4.13) 

Here the quantities F l are functions of £, the quantities listed in (3.16) and (3.17), (3 A 
and their spatial derivatives, t\a and the first derivatives of W and V with respect to t 
and x. They depend linearly on the derivatives of W and V. The mass shell condition 
v° = a/1 + 5ijV l vi defines v° in terms of v\ The characteristics of the equation (4.13) 
satisfy the system: 

dx'/ds = [aA-^v^v ) - p 1 ^ 
dv l /ds = F l 



(4.14) 



The spacetimes considered here have two local Killing vectors. If k a is a Killing vector in 
any spacetime and p a the unit tangent vector to a timelike geodesic, then the quantity p a k a 
is conserved along the geodesic. This allows two conserved quantities for the equations 
(4.14) to be derived. They can be computed, using (4.7) and (4.8) to be: 

Aae w/2 [cosh(V/2)v 2 + sinh(Vy2)i> 3 ] 

Aae- w/2 [sinh(V/2)v 2 + cosh(V72)?; 3 ] ' " '"' 

It is easy to solve for v 2 and v 3 in terms of these two conserved quantities and the bound- 
edness of W and V implies that v 2 and i> 3 are bounded along a characteristic. Consider 
now a solution of the Einstein-Vlasov system where the initial datum for the distribu- 
tion function has compact support. Let P(t) be the supremum of \v\ over the support of 
f(t). Since a, A, W and V have already been controlled pointwise the components Tab 
of the energy-momentum tensor occurring on the right hand side of (4.1) and (4.2) can 
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be estimated in terms of the corresponding frame components. Looking at the explicit 
expressions for these frame components and using the boundedness of v 2 and v 3 in the 
support of / shows that: 

\\T A B{t)\\oo<CP{t) (4.16) 

where C is a constant which only depends on the initial data. To make use of (4.16) an 
estimate for v 1 must be obtained. Define: 



Q(t) = II^W^IU + || W*)||oo + || W*)||oo + IIW*)lloo (4.17) 

Lemma 4.2 If the hypotheses of Theorem 3.1 are satisfied by a solution of the Einstein- 
Vlasov system then the following inequality holds for t^ < t 3 : 

(rtz-ti \ 
1 + P(t 3 ) + J l + P(t 3 -t)+Q(t 3 -t)dtj (4.18) 

The analogous inequality holds for t^ > t 3 . 

Proof In a 3+1 decomposition of a general spacetime the Vlasov equation takes the follow- 
ing form when expressed in terms of frame components: 

df/dt + (avyv°e? - P a )df/dx a 

- [ ei (a)v° + a(-k ab ete b J + 7 £ J V + v k / v°]d f / dv l = 

Here 7^ and 7j fc are Ricci rotation coefficients. Consider the terms appearing in F 1 in 
the case of the symmetry considered here which contain derivatives of V and W. No 
such terms arise from the terms in (4.19) involving the derivatives of a and the second 
fundamental form. To go further it is necessary to have more explicit expressions for 
the rotation coefficients. The four-dimensional rotation coefficients Yjk are identical with 
corresponding three-dimensional ones while: 

7o, = -a-W kj 0aP a + \a-\e^ a (3 h di - S is ey a p%S jt + c) - V s c% t ) (4.20) 

Here = e a ^d t Q l a - Each term in the expressions for the coefficients of the Vlasov equation 
is either independent of the derivatives of W and V, in which case it is bounded as a 
consequence of the estimates already proved, or it is linear in these derivatives. Consider 
now the equation for v 1 in the characteristic system. The estimate (4.18) is obtained by 
considering the dependence on v of those terms which are linear in the derivatives of W 
and V. The quantity v % is bounded unless i = 1 while the quantity v^v k /v is bounded 
unless j = k = 1. However, by the symmetry properties of the rotation coefficients, 
7^ = 0. Hence all terms on the right hand side of the equation for v 1 can be bounded by 
an expression of the form C(l + P + Q). □ 

Lemma 4.3 If the hypotheses of Theorem 3.1 are satisfied by a solution of the Einstein- 
Vlasov system then the quantities P, dtW , dtV d x W , d x V , p, a -1 , the derivative with 
respect to x of all the quantities in (3.16) and (3.17), d x r]A and d 2 (3 A are bounded. 
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Proof The first step is to obtain an estimate for the first derivatives of W and V. In order 
to do this, it is useful to write equations (4.1) and (4.2) in a slightly different way. 



V a V a W + tanhVV a WV a V = -(2/r)V a rV a W - (cosh V)- 1 [e- w T 22 - e w T 33 

- \{e~ w { m f - e w {^f)] (4.21) 
V a V a V - sinhF coshVV a WV a W = -(2/r)V a rV a V - 2(coshy)- 1 [(T 23 - fh AB T AB g 23 ) 

-i(V2V3-i(h AB VAVB)923)} (4.22) 

The advantage of this is that if the right hand sides of (4.21) and (4.22) are replaced by 
zero the resulting equations are those for a wave map (hyperbolic harmonic map) with 
target space R 2 , endowed with the metric cosh 2 VdW 2 + dV 2 . This is a representation of 
the standard metric of the hyperbolic plane in a certain coordinate system. It is natural 
to try to generalize estimates which have been used in the study of wave maps to the 
present situation. Here this is done with an estimate of Gu [11], who used it to prove 
global existence of classical solutions in the Cauchy problem for wave maps defined on 
two-dimensional Minkowski space. Define two null vectors on the two-dimensional space 
coordinatized by t and r by 

e+ = a^id/dt - f3d/dx) + A^d/dx 

(4.23) 

e_ = a - 1 {d/dt-pd/dx)-A- 1 d/dx 
The (2-dimensional) covariant derivatives V e _e_|_ and V e+ e_ are given by: 



V e _e + =a (6 ++ e + + 6 + _e_) 
V e+ e_ = a _1 (6 ye- + b e_) 



(4.24) 



for some bounded functions &++, , b |_ and b The normalization chosen for the 

vectors e + and e_ here is important, since otherwise the covariant derivatives could contain 
the time derivatives of a and (3 1 , quantities which have not yet been shown to be bounded. 
Let E + and E- be the images of e+ and e_ under the wave map, i.e. 

E+ = e + (W)d/dW + e + (V)d/dV 

E_ = e-(W)d/dW + e-(V)d/dV [ ' ' 

Let 71 and 72 be integral curves of e_ and e + respectively and let 7$ = (Wo^, Vo^). The 
observation of Gu is that the equations obtained from (4.21) and (4.22) by replacing the 
right hand side by zero and the given metric by the flat metric just say that E+ is parallelly 
transported along 71 and that E_ is parallelly transported along 72. A similar calculation 
can be done for the equations (4.21) and (4.22) and this gives rise to the following equation 
along 71 (and an analogous equation along 72): 

V aE _E + = (b ++ - r- 1 ae_(r))E + + (&+_ - r~ 1 ae + {r))E_ + B_ (4.26) 
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where £?_ satisfies an inequality of the form \B_\ < C(l + ||Iab||oo)- These equations 
allow the lengths of the vectors E + and E- to be controlled. Multiplying (4.26) and the 
analogous equation for E + by a allows the following inequality to be derived: 

Q(U) < C[Q(t 3 ) + [ 3 4 l + Q(t3-t) + ||T AB (t 3 -t)||oo^] (4.27) 
Jo 

Putting together (4.16), (4.18) and (4.27) gives: 

(1 + P + Q)(U) < (1 + P + Q)(t 3 ) + C [ 3 \l + P + Q)(t 3 -t)dt (4.28) 

Jo 

Hence by Gronwall's lemma P, dtW, dfV, d x W and d x V are bounded. It then follows 
immediately that p is bounded and (3.10) shows that a -1 is bounded. The equations 
(3.9)-(3.14) can be used directly to show that the first derivatives with respect to x of 
all the quantities in (3.16) and (3.17) are bounded. It follows from (4.11) that d x r\A is 
bounded and from (4.12) that d x f3 A is bounded. □ 

Lemma 4.4 If the hypotheses of Theorem 3.1 are satisfied by a solution of the Einstein- 
Vlasov system then the second derivatives of W and V and the first derivatives of f are 
bounded. 

Proof 'If / were zero (the vacuum case) then it would be rather simple to prove this theorem, 
since the equations obtained by differentiating the equations for W and V with respect to 
x are linear in the highest order derivatives in that case. With the coupling to / things are 
less straightforward. When the Vlasov equation is differentiated with respect to x terms 
come up which involve second derivatives of W and V multiplied by first derivatives of /. 
In other words, there are terms which are quadratic in the quantities to be estimated and 
this precludes a direct application of Gronwall's inequality. This problem can be solved 
using a device of Glassey and Strauss [10], which can be seen in a particularly simple form, 
adequate for the present application, in the paper [8] of Glassey and Schaeffer (see also 
[9]). The equation for W can be written in the following form: 

rV a {n b V b W) = {Yi{W, V)l a + Y 2 (W, V)n a )V a W + Z(W, V) (4.29) 

where Z{W,V) contains no derivatives of W or V and Yi(W, V) and Y 2 {W,V) contain 
them at most linearly. Here, for ease of notation, / = e_|_ and n = e_. The equation for V 
can of course be written in a similar form. There are also alternative forms of both of these 
equations where the roles of / and n are interchanged. Differentiating equation (4.29) with 
respect to x gives an equation of the form 

l a V a (d x (n b V b W)) = (Y 1 (W, V)l a + Y 2 (W, V)n a )d x (V a W) + Z(W, V) (4.30) 

where the expression Z(W, V) does not depend on second derivatives of W and V. Suppose 
now that we integrate the equation (4.30) along the characteristic which is an integral curve 
of l a . The only terms which cannot be bounded straightforwardly (even before integration) 
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are those which contain derivatives of the energy momentum tensor with respect to x. It 
will now be shown how a typical term of this type can be handled. The others which occur 
can be taken care of in a strictly analogous way. 
The term which is to be bounded is: 



I" \{e w coshV)- 1 d x T 22 ]{t s -t)dt (4.31) 
Jo 

In fact the coordinate components of the energy-momentum tensor may be replaced by 
frame components at this stage since their spatial derivatives only differ by terms which are 
bounded. Substituting the definition of the frame component T(e 2 , e 2 ) into the expression 
of interest gives 

U J[(e w coshV)- 1 ^) 2 ^ + \v\ 2 )- 1 ]d x fdvdt (4.32) 

The idea of [10] is to express d x as a linear combination of / and the vector 

m = d/dt + (aA-^v^v ) - ^)d/dx (4.33) 

The result is: 

d/dx = a- 1 A(l-v 1 /v°)- 1 (l-m) (4.34) 

This allows the integral in (4.32) to be rewritten as a sum of two terms, one containing I and 
the other containing m. Now it is possible to substitute for mf using the Vlasov equation 
and the result contains only derivatives of / with respect to the velocity variables. These 
derivatives can be eliminated by an integration by parts in v and the result is a bounded 
quantity. The other term is equal to 

t3 t \a- 1 A(e w coshV)-\v 2 )\l + \v\ 2 )- 1 rV a f]( 7l (t 3 -t))dt(l-v 1 /v°)- 1 dv 

= J [(a~ 1 A(e w coskVy^it^x^fit^x^v) 

- (a- 1 A(e w coshF)" 1 )^, x 5 )f(t 3 , x 5 , v)](l - v'/v^dv + ... 

(4-35) 

Thus it is also bounded. The same trick can be applied when W is replaced by V and 
when / and n are interchanged. (In the last case d/dx must be replaced by a combination 
of m and n.) The result of all this is that if 

Qi = HWVaWOHoo + R(n a V a W)|U + ||(/ a V a y)||oo + ||(n°Va^)||oo (4-36) 
then an estimate of the form: 

i + Qi(* 4 ) < i + Qi(h) + C f 3 V + Qife-t))^ (4.37) 
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is obtained. It follows from GronwalPs inequality that Qi is bounded. Hence the derivatives 
W xx , Wtx, V xx and Vt x are bounded. Using this information in the equations obtained by 
differentiating the Vlasov equation with respect to x or v shows that the first derivatives 
of / with respect to these variables are bounded. □ 

5. The main result 

In this section the estimates collected in Section 4 are applied to prove the main result. 
First one last auxiliary lemma is required. 

Lemma 5.1 Consider a CMC initial data set for the Einstein- Vlasov system with local 
{7(1) x i7(l) symmetry. Then there exists a local Cauchy evolution of this data which has 
local {7(1) x {7(1) symmetry, so that the hypotheses of Lemma 2.1 are satisfied. Consider 
next a family of initial data sets of this type on the same manifold such that: 

(i) the data in the family are uniformly bounded in the C°° topology 

(ii) the metrics are uniformly positive definite 

(Hi) the supports of the distribution functions are contained in a common compact set 
(iv) the mean curvatures are uniformly bounded away from zero 

Then the time interval in the conclusion of Lemma 2.1 can be chosen uniformly for the 
Cauchy evolutions of all data in the family. 

Proof The first statement of the proof is essentially a direct consequence of the standard 
local existence theorem for the Einstein- Vlasov system and for CMC hypersurfaces and 
the fact that the resulting spacetimes inherit any symmetry which is present. When 
there is only local symmetry the inheritance of symmetry argument should be applied to 
the universal cover (cf. [17]). The second part of the lemma, concerning families is a 
consequence of the stability of various operations. Firstly, the statement is used that the 
time of existence of a solution of the Einstein- Vlasov system, measured with respect to an 
appropriate time coordinate, depends only on the size of the initial data and that on a fixed 
closed time interval the solution depends continuously on the initial data. Secondly, the 
fact is used that the interval on which a CMC foliation exists in a neighbourhood of a given 
CMC hypersurface depends only on the size of the metric coefficients in an appropriate 
coordinate system and a positive lower bound for the lapse function, provided the mean 
curvature of the starting hypersurface remains bounded away from zero. □ 

Theorem 5.1 Let (M,g,f) be a C°° solution of the Einstein- Vlasov system with local 
{7(1) x {7(1) symmetry which is the maximal globally hyperbolic development of data on 
a hypersurface of constant mean curvature Hq < 0. Then the part of the spacetime to 
the past of the initial hypersurface can be covered by a foliation of CMC hypersurfaces 
with the mean curvature taking all values in the interval (—oo,Hq]. Moreover, the CMC 
foliation can be extended to the future of the initial hypersurface in such a way that the 
mean curvature attains all negative real values. 

Proof Let T be the largest number (possibly infinite) such that the local foliation by CMC 
hypersurfaces which exists near the initial hypersurface can be extended so that the mean 
curvature takes on all values in the interval (—T,Hq). Suppose that T is finite. Then 
Theorem 3.1 and the results of Section 4 imply the boundedness of many quantities on 
the interval (— T, Hq]. It will now be shown by induction that the spatial derivatives of 
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all orders of all quantities of interest are bounded on the given interval, the inductive 
hypothesis is that the following quantities are bounded: 

D n f, D n+1 W, D n (8 t W), D n+1 V, D n (d t V),D n+1 a, 

D n+1 f3\D n+1 A, D n (d t A),D n K, D n rj A , D n+1 f3 A ^'^ 

It follows from the results of Section 4, and in particular Lemma 4.4, that the inductive 
hypothesis is satisfied for n = 1. Suppose now that it is satisfied for a given value of n. 
Then it follows immediately from the field equations (3.9)-(3.14) and (4.11)-(4.12) that 
all quantities which are required to be bounded by the inductive hypothesis at the next 
step are bounded, except possibly for the relevant derivatives of /, W and V. Consider 
the equation obtained by differentiating the Vlasov equation n + 1 times with respect to 
x. There results a linear equation for D n+1 f with coefficients which are known to be 
bounded, except for terms involving derivatives of W and V in the inhomogeneous term. 
If F n {t) HPVIloo and 

Q n (t) = ||£> n+1 ^||oo + ||D n+1 y||oo + P^WIloo + PnWIU 

then this equation implies an inequality of the form 

F n +i (t) < F n+1 (t ) + C / F n+1 (t -s) + Qn+i (t - s)ds (5.2) 

JO 

Similarly, differentiating equations (4.29) and (4.30) n times with respect to x gives a linear 
system of equations for derivatives of V and W with coefficients which are known to be 
bounded, except for terms involving derivatives of order n + 1 of matter quantities in the 
inhomogeneous term. Hence: 

G n+1 (t) < G n+1 (0) + C [ ° F n+1 (t -s) + G n+1 (t - s)ds (5.3) 

Jo 

Putting together (5.2) and (5.3) and applying Gronwall's inequality proves that F n+ i and 
Gn+i are bounded and completes the inductive step. Thus the quantities in (5.1) are 
bounded for all n. Consider now the data obtained by restricting the given solution to the 
hypersurfaces t =const. By what has just been proved, this family of data satisfies the 
conditions of Lemma 5.1. Hence there exists some e > such that each of these initial data 
has a corresponding solution on a time interval of length e. Hence the original solution 
extends to the interval (— T — e,Ho), contradicting the maximality of T. It follows that 
in fact T = oo, as desired. This means in particular that the spacetime has a crushing 
singularity in the past, and hence that the CMC foliation covers the entire past of the 
initial hypersurface. 

Now let T' be the largest number such that the CMC foliation can be extended to the 
interval (— oo, T"). Since the spacetime contains no compact maximal hypersurface T' < 0. 
If T' were strictly less than zero it could be argued as in the first part of the proof that 
the CMC foliation could be extended further, which would contradict the definition of T' . 
Hence in fact T" = 0. □ 
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This argument does not prove that the entire future of the initial hypersurface is covered 
by the CMC foliation. In connection with this it is interesting to note that if instead of 
assuming, as is done in this paper, that the cosmological constant A vanishes, it is assumed 
that A < then the same types of arguments apply to give a stronger theorem. (The choice 
of sign convention for the cosmological constant used here is such that A < corresponds 
to anti-de Sitter space.) With A < the result is that the whole spacetime can be covered 
by a CMC foliation with the mean curvature taking on all real values. The reason for this 
difference can be traced to the estimate for a following from the lapse equation, which in 
general reads a < (\t 2 — A) -1 . 

6. The case of wave maps 

In this section we consider what happens when the collisionless matter described by 
the Vlasov equation is replaced by a wave map as source in the Einstein equations. This is 
quite natural, given that, as was seen in Section 4, a wave map comes up automatically in 
the case of vacuum spacetimes. Let (N, h) be a complete Riemannian manifold. If (M, g) 
is a Lorentz manifold a wave map from M to iV is a map which satisfies the equation 
whose expression in local coordinates x a on M and y 1 on N is: 

v a vV + r^v Q </> J v a <^ = o (6.i) 

(Wave maps are also known as (hyperbolic) harmonic maps or nonlinear sigma models.) 
The global Cauchy problem for wave maps on two-dimensional Minkowski space was solved 
by Gu[ll] and for wave maps on three-dimensional Minkowski space which are invariant 
or equivariant under rotations by Christodoulou, Shatah and Tahvildar-Zadeh [3, 4, 19]. 
The results of [4] were applied to the Einstein-Maxwell equations in [2]. Associated to a 
wave map <f> is the energy- momentum tensor: 

T a(3 = [V a 7 V^0 J - ±(V 7 7 VV)<7a/3]/*/J (6.2) 

and this can be used to couple the wave map to the Einstein equations. In harmonic 
coordinates the coupled equations form a system of nonlinear wave equations and so a 
local existence and uniqueness theorem can be proved by the usual methods. The energy- 
momentum tensor of a wave map satisfies both the dominant and strong energy conditions. 
This can be seen by noting that both these conditions are purely algebraic in nature and 
can be checked using normal coordinates based at a given point of N. Then the energy- 
momentum is reduced at a point to a sum of terms, each of which is the energy-momentum 
tensor of a massless scalar field. 

Consider now the case of a solution of the Einstein equations with local U(l) x U(l) 
symmetry coupled to an invariant wave map. To say that the wave map is invariant means 
that each surface of symmetry is mapped to a single point of N. Since the relevant energy 
conditions hold, it follows that the analogues of the results obtained in Sections 2 and 3 
for the Einstein-Vlasov system are also valid for the Einstein-wave map system. Given 
that in proving Theorem 5.1 a wave map was already estimated, albeit for a special target 
manifold (AT, h) , it appears straightforward to generalize that theorem to the case of the 
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Einstein-wave map system. In fact the equation of motion for the wave map does not 
involve W, V or t\a while the combinations of matter terms occurring in (4.1) and (4.2) 
vanish identically for the energy-momentum tensor of an invariant wave map. Thus there 
is no direct coupling between the wave map describing the matter and the wave-map-like 
equation satisfied by W and V. The one difficulty which occurs is that, in contrast to 
the special case of the hyperbolic plane, there is no global coordinate system on N in the 
general case. The equation for an invariant wave map can be written in the form: 

V a (r 2 V V) + r 2 r^V a J V o 0* = (6.3) 

This bears a strong resemblance to equations (4.1)-(4.2), with the difference that there 
are no terms involving rj or matter quantities in (6.3). This makes the analogue of the 
calculation (4.9) for the wave map superfluous. This is just as well, since it seems difficult 
to formulate an analogue of (4.9) in the case that there is no global coordinate system on 
N. What can be done instead is to go directly to the analogue of (4.26) for the wave map. 
Define: 

E+ = e+U^d/dcf) 1 

~ T (6-4) 

Then E + and E- satisfy propagation equations like (4.26) along 71 and 72 respectively. 
There is no term corresponding to B- in this case. It follows that under the hypotheses of 
Theorem 3.1 the length of the vectors E+ and E- is bounded on the given time interval. 
This implies a bound on the distance of any point of the image under of this time interval 
from the image of the initial hypersurface. In particular, the image of this time interval 
under is contained in a compact subset of N. This compact set can be covered by a finite 
number of charts, each of which can be chosen to be defined on a domain with compact 
closure in a larger chart domain. In each of these charts the quantities </> J , dt4> T and <9 X </> J 
are bounded. Moreover, any of these charts the following analogue of (4.30) holds: 

Z o V a (0*(ra 6 V 6 J )) = (Y 1 (<p I )l a + Y 2 { ( j> I )n a )d x {V a (j> 1 ) + Z^ 1 ) (6.5) 

This equation and the equations obtained by differentiating it repeatedly with respect to 
x can be used to inductively bound all spatial derivatives of (j) 1 . This proceeds essentially 
as in the proof of Theorem 5.1; it is merely necessary to be careful about the different 
charts which occur. In the case of a wave map define F n (t) to be the maximum over 
the finite set of charts chosen of ll-D"" 1 " 1 ^!^ and ||-D n (<9t0 7 )||oo- When the derivatives of 
lower orders are known to be bounded, this is equivalent to choosing for each point one 
chart which contains its image and only taking the supremum over those values. When 
the quantity F n (t) is bounded the derivatives of order n of the frame components of the 
energy-momentum tensor are bounded. In order to get an inequality which can be used to 
control F n (t), we would like to integrate a derivative of (6.5) along a characteristic (integral 
curve of lor n). The image of this characteristic under <p need not be contained in a single 
chart. Consider such a characteristic 7, parametrized by t from t = to t = T. For each 
t E [0, T] there exists an interval /, open in [0, T], whose image under 4> o 7 is contained in 
one of the chosen charts on N. By compactness of [0, T], finitely many of these intervals 
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cover it. It follows that there is a finite sequence of times {0 = t±, ti, ■ ■ ■ , tk} = T such 
that 4> o 7([tj, U+i]) is contained in one of the chosen charts for all i between 1 and k — 1. 
It follows from (6.5) that 

F n (t k ) < Fnitk-^e ^"-^ (6.6) 

This is enough to allow F n to be bounded for all t G [0, T]. Thus the following analogue 
of Theorem 5.1 is obtained: 

Theorem 6.1 Let (N,h) be a complete Riemannian manifold and let (M,g,(j>) be a C°° 
solution of the Einstein equations with local U(l) x U(l) symmetry coupled to an invariant 
wave map with target space (N, h) which is the maximal globally hyperbolic development of 
data on a hypersurface of constant mean curvature H < 0. Then the part of the spacetime 
to the past of the initial hypersurface can be covered by a foliation of CMC hypersurfaces 
with the mean curvature taking all values in the interval (—oo,Hq\. Moreover, the CMC 
foliation can be extended to the future of the initial hypersurface in such a way that the 
mean curvature attains all negative real values. 
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